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Abstract 

Using a large- expansion we present and solve the linearized equation for 
the superconducting gap for a generalized t — J model which also contains 
phonons within a Holstein model. Keeping all terms up to 0{1/N) the kernel 
of the gap equation consists of an electron-phonon part with self-energies and 
vertex functions renormalized by the interactions of the t — J model, and a 
t — J part unaffected by phonons. Considering first the electron-phonon part 
we find that the leading superconducting instability always occurs in the s- 
wave channel with a Tc which is lower compared to that of non-interacting 
electrons, especially at large dopings. The corresponding isotope coefficient 
(3 is larger than 1/2 and increases with decreasing doping. Including also the 
t — J part in the gap equation the leading Tc has always d-wave symmetry with 
phonons giving a positive contribution to Tc. [3 is very small at small dopings 
and increases towards the classical value 1/2 with increasing dopings similar 
as in many cuprates. Considered as a function of the phonon frequency ojq at 
a fixed coupling strength f3 increases monotonously from zero to about 1/2 
with increasing u)q. 

PACS numbers: 74.20.Mn, 74.25.Kc, 74.20-z 
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I. INTRODUCTION 



High-Tc cuprates show strong electronic correlations but also a non-neglegible electron- 
phonon interaction. The latter causes, for instance, the observed isotope effect in the transi- 
tion temperature T^, especially, away from optimal dopingB. It is thus of interest to calculate 
superconducting instabilities taking both the electron-electron and the electron-phonon in- 
teraction into account. There are several calculations dealing with the case of weak electronic 
correlationsoQ. In these calculations the phonon-mediated part of the effective interaction 
is unaffected by electronic correlations. The electronic part to the effective interaction is 
calculated using a Hubbard model and assuming U to be small compared to the band width. 
It is repulsive and strongly peaked near the M-point in the Brillouin zone. As an alternative, 
Refs.ili use the spin-fluctuation term of Ref.i. Since the phonon part is attractive through- 
out the Brillouin zone it is plausible that the d-wave pairing due to the electronic part will 
be weakened if the electrons couple strongly to phonons near the M— point and only weakly 
affected if only phonons with small momentum are involvedlli^'i. Refs.i'@, for instance, find 
that the for d-wave superconductivity is always lowered by phonons and that the isotope 
coefficient f3 is negative for a constant electron-phonon coupling function. If only phonons 
with small momenta play a role (3 becomes positive!. Ref.i even finds a change of sign in (3 
as a function of the phonon frequency. 

The case of strong electronic correlations has quite different features compared to the 
weak correlated case. The purely electronic part of the effective interaction is rather a smooth 
function of momentum without a sharp peak at the M-point due to spin fluctuations0. 
Furthermore, electronic correlations modify substantially the phonon-mediated part of the 
effective interaction0Hll. It is the purpose of this communication to present results for this 
case of strong electronic correlations. A suitable model for such calculations is a generalized 
t — J model which also includes phonons within a Holstein model. Extending the spin 
degrees from 2 to systematic approximations for the effective interaction in terms of 
powers in 1/A^ can be carried outJiirtJ. In a first step we consider only the phonon-mediated 
effective interaction which, unlike in the case of weak correlations, is strongly modified 
by vertex corrections transforming the original constant electron-phonon coupling of the 
Holstein model into a strongly momentum- and frequency-dependent function. Explicit 
results for Tc and (3 will be given for the leading symmetry channels of the superconducting 
order parameter as a function of the doping. In a second step these results are extended to 
the case where also the purely electronic contribution to the effective interaction is taken 
into account. 



II. LINEARIZED EQUATION FOR THE SUPERCONDUCTIVITY GAP 

Our Hamiltonian for the t — J model plus phonons can be written as 

AT t J A AT « J 

P — 1...N p,q — l...N 



. 4Ar 

i,q = l...iV 
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Y: -^[al + a,](X,^^-<X,^^>). (1) 



The first three terms correspond for = 2 to the usual t — J modeL For = 2 X is 
identical with Hubbard's projection operator Xf = |^ >< '^|, where |^ > denotes for p = 
an empty and for p = 1,2 a singly occupied state with spin up and down, respectively, at the 
site i. tij and Jij are hopping and Heisenberg interaction matrix elements between the sites 
The fourth term in Eq.(l) represents one branch of dispersionless, harmonic phonons 
with frequency Uq. The fifth term in Eq.(l) describes a local coupling between the phonon 
and the change in the electronic density at site i with the coupling constant g. < X > 
denotes the expectation value of X. The extension from N = 2 in Eq.(l) to a general N 
has been discussed in detail in Ref.0. The label p runs then not only over the two spin 
directions but also over N/2 identical copies of the orbital. The symmetry group of H is the 
symplectic group Sp{N/2) which allows to perform 1/N expansions for physical observables. 
In Eq.(l) the electron-phonon coupling is scaled as 1/Vn whereas the free phonon part is 
independent of N. As a result the leading contributions to superconductivity from the t — J 
model alone and from the phonons are both of order 0{1/N) which allows to treat them on 
an equal footing. 

Instabilities towards superconductivity in the t — J model have been studied in the 
above framework in Refs.lljt3il£l. The contributions to the anomalous self-energy from the 
phonon-mediated effective interaction have been derived in Refs.0'0 for the case Jij = 0. 
Generalizing the latter treatment to a finite J similar as in Refill the linearized equation 
for the superconducting gap E^^ for the entire Hamiltonian Eq.(l) can be written as 

T^anik) = — y e(k, k')—^ \-—T.an{k'). (2) 

Nc is the number of cells and k the supervector k = (u;„,k), where u;„ denotes a fermionic 
Matsubara frequency and k the wave vector. e(k) is the one-particle energy which is un- 
changed by the phonons in 0(1) and thus given by Eq.(41) of Ref.0. The kernel 6 in Eq.(2) 
consists of two parts 

e(k, k') = e'-\k, k') + Q'-'p{k, k'). (3) 

The first contribution G*~'^ in Eq.(3) comes from the t — J model. Explicit expressions for 
it have been given in Eqs.(42)-(52) in Ref.0. The second term 9*^ ^ in Eq.(3) is due to 
phonon-mediated interactions and is given by 

e-^(A;, k') = ^^'""^ 



■ ^,{k',k-k'h,{k,k'-k). (4) 

7c is the charge vertex for which an explicit expression has been derived in Ref.^. For uncor- 
related electrons there are no vertex corrections, i.e., 7c = —1. Since g and ujq are assumed 
to be independent of k 6"^"^ is also independent of k and only s-wave superconductivity is 
possible. When correlations are present the kernel Q'^'P depends on k and Un through the 
charge vertex 7c and general symmetries for the order parameter may become possible. 
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III. CALCULATION OF Tc AND (3 FROM THE PHONON-MEDIATED PART 



Taking a square lattice with point group C^^ Q{k,k') and e(k) in Eq.(2) are invariant 
under C^^ which means that can be classified according to the five irreducible repre- 

sentations Fj of C^y. s-wave symmetry corresponds to Fi, d-wave symmetry to F3, etc. In 
the weak-coupling case Q{k,k') can be approximated by its static limit 0(k, k'). Putting 
all momenta right onto the Fermi line the sum over k' in Eq.(2) can be transformed into a 
line integral along the Fermi line. Assuming a certain irreducible representation Fj for the 
order parameter, the hne integral can be restricted to the irreducible Brillouin zone (IBZ) 
introducing a symmetry-projected kernel ©^(k, k') with k, k' e IBZ. Finally, the hne integral 
was discrctizcd by a set of points [k^] along the Fermi line in the IBZ with line elements 
[s(k^)]. Denoting the smallest eigenvalue of the symmetric matrix 



.(k^).(k 



:FV\®^{K.^) (5) 



4^^^|V6(k^)|-|V6(k^, 

by Aj Eq.(2) yields for \ N — 2, and in the weak-couphng case the BCS-formula 

= 1.13a;oe'/^^ (6) 

As usual we took uoq as a suitable cutoff. If Aj > we have, of course, = 0. According to 
Eq.(6) the absohitc vahic of Aj characterizes the strength of the effective interaction in the 
symmetry channel i. The overall strength of the electron-phonon coupling is conventionally 
expressed in terms of the dimensionless coupling A defined by 

A = — • (7) 

8a;o ^ ' 

In Eq.(7) we have introduced a factor 1/2 to account for the prefactor 1/2 in Eq.(2) after 
setting there N = 2. We also used the average density of states l/(8t) for the density of 
states factor in A and put t equal to 1. For the range of dopings we will be interested in, 
i.e., 0.15 < 6 < 0.8, the density of states varies only little with doping so that the definition 
Eq.(7) of A is appropriate for all dopings. 
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FIG. 1. Phonon contribution to the lowest eigenvalues Aj of the static kernel @ for the 
representations Fi and F3 of 6*4^ as a function of the doping 6. 

Fig. 1 shows the eigenvalues Aj for the two leading symmetries i = 1,3 using the 
parameter value J = 0.3 and t as the energy unit. 6 is the doping measured from half-filling. 
For A we have chosen the value 0.75 which is a typical value obtained in LDA calculations0'0. 
The lowest eigenvalue occurs always in s-wave symmetry, i.e., for i = 1. It depends only 
weakly on doping for 0.3<5<0.8. In this region the order parameter also varies only 
slowly along the Fermi line, describing thus rather isotropic s-wave superconductivity. Below 

6 = 0.3 Ai rapidly decreases with decreasing doping wich is caused by a soft mode which 
freezes into an incommensurate bond-order wave at S bo ~ 0.140. As a result the s-wave 
order parameter becomes less and less isotropic with decreasing 6 changing, for instance, 
for 6 = 0.2 by about a factor 3 along the Fermi line but it does not pass through zero. 
The d-wave coupling constant A3 would be zero in the uncorrelated case. Decreasing 6 from 
large values correlation effects increase and the static vertex function 7(k, k — k'ldevelops 
more and more a forward scattering peak in the transferred momentum k — k'0@. In the 
extreme case where 7 is proportional to S(k — k') Ai and A3 would become degenerate. Fig. 
1 shows that this degeneracy is nearly reached at low dopings. The strong localization of 
the effective interaction is caused by the forward scattering peak in 7 and, in addition, by 
the incipient instability at 6bo which is also very localized in k-space. 

The neglect of retardation effects which leads to the BCS-formula Eq.(6) is doubtful 
for several reasons. First, A = 0.75 no longer corresponds really to a weak-coupling case. 
Secondly, is no longer small compared to the frequency of the soft mode near 6bo which 
means that the frequency dependence of cannot be neglected for these dopings. Thirdly, 
and most severely, the formation of the forward scattering peak in the static vertex function 

7 is accompagnied by a strong frequency dependence which should be taken into account 
on the same footing as its momentum dependence. We thus have solved the gap equation 
Eq.(2) assuming only that the momenta can be put right onto the Fermi line but keeping 
the full frequency and momentum dependence along the Fermi line. 
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FIG. 2. Transition temperature Tc in units of t as a function of the doping 6 for the uncorrelated 
(dash-dotted hne, Fi symmetry) and correlated (dashed and sohd hues for Fi and F3 symmetries, 
resp.) cases, using only the phonon contribution. 

Fig. 2 shows the obtained results for for a s-wave (dashed line) and a d-wave (solid 
line) order parameter, together with the curve for the uncorrelated case (dash-dotted line). 
The first observation is that correlation effects always suppress s-wave superconductivity. 
The suppression is large for 6 > 0.25 and becomes less pronounced at smaller dopings in 
agreement with previous results based on the static approximation^. The order parameter 
associated with the dashed line in Fig. 2 varies only slowly along the Fermi line, i.e., 
we have a usual s-wave order parameter without nodes. In a Gutzwiller description this 
order parameter would be exactly zero if retardation effects can be neglected and N = 2 
is taken: Integrating out the phonons the effective interaction becomes in the static limit 
proportional to the double occupancy operator. Any matrix element formed with Gutzwiller 
wave functions thus would be zero. In contrast to that the enforcement of our constraint at 
large A^'s, namely, that only A^/2 out of the total states at a given site can be occupied 
at the same time, gives rise to another effect, which is absent in the Gutzwiller treatment: 
the effective interaction becomes more and more long-ranged with decreasing doping. This 
makes isotropic s-wave superconductivity possible even in the presence of a A^ = 2 constraint. 
Fig. 2 nevertheless shows that the reduction of due to the constraint is substantial except 
at very small dopings where the effective interaction becomes extremely long-ranged and 
the suppression by the constraint is small. 

The kernel B*^"^ is in the absence of correlations and in the static limit negative for 
all arguments k and k'. This implies that the s-wave symmetry has always the highest Tc. 
According to our numerical studies the same is also true in the correlated case which explains 
why the solid line in Fig. 2, describing d-wave superconductivity, is always below the dashed 
line. As indicated previously the very existence of finite values for with d-wave symmetry 
is in our model due to electronic correlation effects. Fig. 1 showed that the momentum 
dependence of 7c causes in the static limit a finite coupling strength in the d-wave channel. 
The solid line in Fig. 2 proves that a finite Tc in the d-wave channel results from this even 
if the strong frequency dependence of 7c is also taken into account. 
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FIG. 3. Isotope coefficient /3 as a function of the doping 5 considering only the phonon 
contribution for Fi (dashed line) and F3 (solid line) symmetries. The filled circles are calculated 
values. 

Fig. 3 shows results for the isotope coefficient f3 in the case of s-wave (dashed line) 
and d-wave (solid line) superconductivty. We assumed hereby that Uq is inversely and g 
directly proportional to the the square root of the mass rendering A independent of the 
mass. Neglecting the frequency dependence of 7c there is only one energy scale involved in 
the gap equation which yields /3 = 1/2 independent of doping, A, etc. 

In contrast to that the curves in Fig. 3 show that (3 increases monotonically with de- 
creasing doping and is always larger than 1/2. The deviaton of (3 from the value 1/2 is 
due to the frequency dependence of 7c which increases with decreasing doping0 and acts as 
a second energy scale. The non-adiabatic effects produced by electronic correlations thus 
always increase f3 for our range of parameter values. This should be contrasted to the case 
where vertex corrections due to the electron-phonon interaction has been studiedi'§, and 
where f3 may be larger or smaller than 1/2. We can conclude from our results that models 
based on phonon-induced effective interactions in the presence of correlations seem not to 
be very attractive models for high-Tc superconductivity because a) the largest transition 
temperatures are found for s-wave symmetry for all dopings and b) the isotope coefficient is 
larger than 0.5, especially at low dopings. Experimentally, high-Tc cuprates exhibit d-wave 
superconductivity and, at least near optimal dopings, small isotope effects. 




IV. CALCULATION OF Tq AND /3 FROM THE FULL EFFECTIVE 

INTERACTION 

Superconducting instabilites of the pure t — J model have been studied in0§0. The 
smallest eigenvalue of B*~'^ and, correspondingly, the largest Tc always occur in the d-wave 
symmetry channel. On the other hand, the phonon-induced effective interaction described 
by B'^~^ favors in the t — J model mainly s-wave superconductivity. In this section, we will 
consider the case where both effective interactions are present and study the symmetry of 
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the resulting superconducting state, its transition temperature and isotope coefficient. 
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FIG. 4. Lowest eigenvalues Aj of the total static kernel G for the representations Fi and F3 of 
as a function of the doping 6. 

Fig. 4 shows the two lowest eigenvalues of the matrix Eq.(5) which occur in the s- and 
d-wave symmetry channel. For B both terms in Eq.(3) have been included. For 6 < 0.6 the 
lowest eigenvalue has d-wave symmetry and decreases steadily and strongly with decreasing 
doping diverging finally at 6 bo- Comparing the solid lines in Figs. 1 and 4 of this work 
and Fig. 2 intd one sees that the lowest eigenvalue is roughly additive in the t — J and 
the phonon contributions though, of course, B*"*^ and B^~^ do not commute with each 
other. Note also the considerable contribution of the phonon part to A3 though the overall 
dimensionless coupling constant has the rather moderate value A = 0.75. The eigenvalue 
Ai (dashed line in Fig. 4) is rather constant for 6 > 0.3 and decreases much less towards 
lower dopings than A3. There is a crossover from d-wave to s-wave symmetry in the lowest 
eigenvalue at 5 ~ 0.6 suggesting that the s-wave order parameter is more stable than the 
d-wave order parameter at high dopings. Such a crossover is expected for general reasons: 
For large dopings correlation effects play a minor role. As a result the d-wave part in the 
phonon-induced effective interaction as well as the whole t — J part becomes small whereas 
the s-wave part of the Holstein model will dominate. A comparison of the solid lines in Figs. 
1 and 4 nevertheless shows that the latter is very effectively suppressed by the t — J part at 
practically all considered dopings. For instance, |Ai| in Fig. 1 is for 6 > 0.25 about a factor 
4 larger than in Fig. 4. 
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FIG. 5. Transition temperature Tc in units of t as a function of tlie doping 5 for symmetry 
for A = (dasli-dotted line) and A = 0.75 (solid line). 

The solution of the gap equation Eq.(2) is not straightforward because 9 contains both 
instantaneous and retarded contributions and because the various retarded contributions 
have different effective cutoff energies. We developed i a method to solve Eq.(2) directly 
avoiding the use of pseudopotentials. The only simplification is that we put the momenta 
in retarded (but not in instantaneous) contributions to the Fermi line. This approximation 
has been checked numerically and found to be very well satisfied in our case. Using the fact 
that the instantaneous kernel consists only of a few separable contributions Eq.(2) can for a 
given Matsubara frequency reduced to a linear matrix problem of order 1 and 3 for d- and 
s-wave symmetry, respectively. 

In Fig. 5 we have plotted the calculated values for Tc for d-wave symmetry Fa as a 
function of the doping 6. The dashed-dotted line corresponds to A = 0, the solid line to A = 
0.75. The Figure shows that phonons always increase the for d-wave superconductivity. 
The increase is especially large at low dopings. This is quite in contrast to calculations in 
the weak-coupling caseBi where the of d-wave superconductivity is lowered by phonons. 
The physical picture in the latter case is that the effective interaction of the electronic part 
is repulsive and strongly peaked near the M-point whereas the phononic d-wave part is 
attractive and diminuishes especially at the M-point the repulsion causing a lowering of T^. 
In our case the dependence of the static effective interaction of the pure t — J model on the 
transferred momentum can be inferred from Figs. la) and lb) in Ref.ty. For small dopings the 
phonon part is restricted to small momenta which means that it would contribute only near 
the X-point in Fig. lb) of Ref.§. As a result the total effective interaction would become 
even more attractive around the X-point and elsewhere not be changed. This clearly would 
enhance the d-wave part of the effective interaction which is also in agreement with Fig. 
4. From Figs. 1, 4 and 5 it is evident that the large lowering of the eigenvalue A3 due to 
phonons corresponds only to a rather moderate increase in Tc. This means that the use of 
a BCS-formula with a fixed effective cutoff would grossly overestimate the increase in Tc for 
d-wave superconductivity due to phonons. What has been overlooked in such an approach 
is that we deal with at least three energy scales, namely t, J, and ujq. For instance, the 
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instantaneous contribution of the t — J part is characterized by the energy scale t whereas 
the phononic one by ujq. Since t >> cuq it is evident that the phononic part in A3 will 
contribute to much less than the t — J part. 

We were unable to find any finite transition temperatures for i 7^ 3, i.e., for symmetries 
different from the d-wave symmetry.. Taking the accuracy of our calculation into account 
this means that < 0.002 for z 7^ 3. As shown in Fig. 2 the phonon-induced effective 
interactions leads to a considerable Tc for s-wave superconductivity. The t — J part to the 
effective interaction, however, is very repulsive in the s-wave channel prohibiting s-wave 
superconductivity. We cannot exclude that the crossover at 5 ~ 0.6 in Fig. 4 in the lowest 
eigenvalue from d-wave to s-wave could stabilize a s-wave order parameter at large dopings. 
Another possibility is that other symmetries than d- or s-wave become stable at large dopings 
in view of the approximate degeneracy of the eigenvalues of practically all symmetries in the 
pure t — J mo del§. In any case, the corresponding T^'s would be smaller than ~ 0.002 for 
all dopings and thus be rather irrelevant. 

The behavior of the isotope coefficient /3 as a function of doping is another interesting 
test for any theory of high-Tc superconductivity. Fig. 6 shows the calculated (3 for A = 0.75, 
J = 0.3, uJo = 0.06 in the case of symmetry. It is always positive, starting from small 
values at small dopings. With increasing doping it increases monotonously approaching 
values near 1/2 at large dopings. Our calculated curve shows the typical behavior of the 
experimetal data, namely, a small value for /5 at optimal doping and substantial values far 
away from optimal doping]. The curve in Fig. 6 is the result of several 
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FIG. 6. Isotope coefficient /3 as a function of the doping 6 for d-wave-like T3 symmetry. The 
filled circles are calculated values. 

competing effects. We saw in Fig. 3 that correlation effects in the phonon-induced part 
cause a doping dependence of /3 which is opposite to that in Fig. 6, namely, a monotonously 
increasing P with decreasing doping starting at the classical value of 1/2 at large dopings. 
Including also the t — J part the large values of (3 at low dopings are nearly completely 
and at larger dopings partially quenched. Part of this effect may be understood from Fig. 
5. It shows that the relative importance of the phonon contribution increases steadily with 
increasing 5 except at very small dopings. In a simple picture one thus may argue that at 
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small dopings Tc is mainly due to the t — J part. This part decays faster than the phononic 
part with doping so that at larger dopings is mainly due to the phonons. A more realistic 
interpretation of Fig. 6 should, however, also take into account the complex competition 
between the three energy scales t, J, and ujq, pair breaking effects, etc. 

In order to get more insight into the behavior of (3 we have calculated Tc and j3 in the Fs 
symmetry channel for the fixed doping value 5 = 0.20 as a function of the phonon frequency 
cuq. Fig. 1^ shows the result for J = 0.3 and A = 0.75. The solid line in Fig. [^represents the 
value for T^, multiplied with 10. Since we keep A fixed the coupling constant g approaches 
zero in the adiabatic limit Uq ^ which means that in this limit is solely due to the t — J 
part. With increasing phonon frequency T^. 
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FIG. 7. Transition temperature multiplied by 10, and isotope coefficient /? as a function 
of the phonon frequency loq for d-wave-like symmetry. 

increases monotonously, passing over from an initial sublinear to a linear behavior. The 
solid line shows that phonons always increase for d-wave superconductivity. The absolute 
values in the Figure are quite remarkable: Though the employed coupling constant for A is 
rather moderate and though the phonon contribution to in the d-wave channel is entirely 
due to correlation effects Tc increases by nearly a factor 3 for tuo = 0.2 corresponding to about 
the largest phonons in the cuprates. The dashed line in Fig. ^ represents the dependence 
of P on uq. For small ujq (3 is practically zero. (3 increases monotonously with increasing u:q 
and tends to the classical value 1/2 at large phonon frequencies. The dashed line in Fig. 
may be understood in physical terms roughly in the following way. In our model the static 
effective coupling is independent of the phonon mass, i.e., of ujq. [3 is thus determined by an 
effective energy cutoff. For small a;o's this cutoff is mainly given by electronic parameters 
in the t — J part leading to a small value for f3. At large c^o's the phonon contribution to 
Tc, which is only due to electronic correlations, determines mainly the total effective energy 
cutoff causing a value of /? near 1/2. 

In conclusion, we have treated the electron-electron and the electron-phonon interactions 
in a generalized t — J model by means of a systematic 1/A^ expansion and have solved the 
resulting linearized equation for the superconducting gap by reliable numerical methods. We 
found that electronic correlations affect phonon-induced superconductivity in several ways: 
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Instabilities towards d-wave or other symmetries different from s-wave become possible. 
The corresponding transition temperatures are, however, always smaller than that of s-wave 
superconductivity. Moreover, the Tc for s-wave superconductivity is at larger dopings heavily 
and at small dopings somewhat suppressed by correlations. The isotope coefficient (3 is 1/2 
at large dopings but increases with increasing correlations, i.e., decreasing dopings, both in 
the s- and d-wave channel. Including also the t — J part in the effective interaction we found 
that within our numerical accurracy only d-wave superconductivity is stable for dopings 
0.15 < 6 < 0.8 and that the phononic part always increases Tc for the above dopings. The 
t — J part in the effective interaction changes the dependence of jS on doping: f3 assumes now 
small values at low dopings and increases monotonously with doping towards the classical 
value 1/2 at large dopings. Keeping A fixed and varying the phonon mass, i.e., ujq we find 
that Tc and (3 increase monotonously with uq and that /3 varies between zero at small and 
roughly 1/2 at large phonon frequencies within the interval < uq < 0.2. 

Our findings differ in many aspects from the corresponding results based on weak- 
coupling calculationdi^. In these calculations for d-wave superconductivity induced 
either by a f/ or a spin-fluctuation term is always suppressed by phonons. For instance, 
in the fluctuation-exchange approximation Tc drops to zero if the phonon-mediated on-site 
attraction Up becomes comparable to the Hubbard term iM. We treated the opposite case 
where U >> Up and found a different behavior, namely, that phonons enhance the Tc for 
d-wave superconductivity. Our different result is mainly caused by corrections to the bare 
electron-phonon vertex due to the strong electron-electron interaction. This vertex devel- 
ops for not too large dopings a forward scattering peak so that only phonon with small 
momenta can couple to the electrons. As a result the electron-phonon and the t — J contri- 
butions to the gap decouple in k space. Our calculations show that the two contributions 
no longer cancel each other to a large extent but, on the contrary, enhance each other. For 
a momentum-independent bare electron-phonon coupling weak-coupling calculations yield 
a small, often negative value for the isotope coefficient /3 which is rather independent of 
doping. In our case the strong momentum-dependence of the effective electron-phonon cou- 
pling, induced by electronic correlations, causes a strong dependence of /5 on doping: being 
always positive, /3 is small at optimal doping and assumes values of roughly 1/2 at large 
dopings in agreement with the basic features of the experimental data in the cuprates. 

The presented results are accurate at large A^'s because we have taken only the leading 
terms of a 1/N expansion. Phonon renormalizations and vertex corrections due to the 
electron-phonon interaction are of 0(1/A^) and thus have been omitted. On the other hand 
it is known that in the physical case N = 2 anharmonic effects in the atomic potentials and 
the formation of polarons occur if A is about 1 or larger0. This suggests that keeping only 
the leading order of the 1/A^ expansion cannot describe adequately the case = 2 at large 
coupling strengths or if the Migdal ratio uo/t is no longer small. Correspondingly, we have 
shown numerical results for rather moderate values for A and small Migdal ratios. 
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